Fermi arcs and the hidden zeros of the Green's function in the pseudogap state 
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We investigate the low energy properties of a correlated metal in the proximity of a Mott insu- 
lator within the Hubbard model in two dimensions. We introduce a new version of the Cellular 
Dynamical Mean Field Theory using cumulants as the basic irreducible objects. These are used 
for re-constructing the lattice quantities from their cluster counterparts. The zero temperature one 
particle Green's function is characterized by the appearance of lines of zeros, in addition to a Fermi 
surface which changes topology as a function of doping. We show that these features are intimately 
connected to the opening of a pseudogap in the one particle spectrum and provide a simple picture 
for the appearance of Fermi arcs. 
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The origin of the pseudogap persists as one of the lead- 
ing unresolved problems in the physics of the copper- 
oxide high temperature superconductors Q. Since a lot of 
the physics of these systems arise from short range corre- 
lations, cluster extensions of single site Dynamical Mean 
Field Theory (DMFT)0 are ideally suited for this prob- 
lem. In fact, using such methodologies, several groups 
have found [3 that a pseudogap, as evidenced by a sup- 
pression of the density of states at the Fermi level, ap- 
pears near the doped Mott insulator as described by the 
Hubbard model. This effect is caused solely by short- 
range correlations and no long range order or preformed 
pairs need to be invoked. In the present work we present 
an extension of the cluster methodology that allows us to 
identifying the emergence of lines of zeros of the Green's 
function at zero energy (i.e. the Luttinger surface) in ad- 
dition to the quasiparticle poles (i.e. the Fermi surface), 
in the proximity to the Mott transition. These results 
are similar to those found in quasi one-dimensional sys- 
tems by Essler and Tsvelik|6j. The appearance and the 
evolution of a pseudogap in the particle spectral function 
is governed by the topology of these lines. At small hole 
doping, the Fermi surface, i.e. the line of poles, is a hole 
pocket having a Luttinger surface in close proximity. The 
quasiparticle weights along a portion of the Fermi contour 
are suppressed by the proximity of the zero line, generat- 
ing the Fermi arc behavior of the spectral function which 
was identified experimentally [l^ . 

We formulate a new cluster approach based on a re- 
summation of a strong coupling expansion around the 
atomic limit, which generalizes the CDMFT approach 5j. 
We use the notations of Reffl for a general lattice Hamil- 
tonian, 

H = Ho + H, (1) 

i a i^ij a, /3, a', /3' 

where the local, Hq, and the non-local, Hi, terms are 
expressed in terms of Hubbard operators, X"^ 3 . Here 



a, (3, a' and f3' represent single-site states. All the on-site 
contributions, like, for example, the Hubbard U interac- 
tion, are included in X a , while the non-local coupling con- 
stants E?? a 13 can be understood as generalized hopping 
matrix elements and may include hopping terms (iy), 
spin-spin interactions (J%j), or non-local Coulomb inter- 
actions (Vij). 

A cluster DMFT scheme^ maps the lattice model onto 
an effective impurity problem on a real space cluster C, 
defined by the statistical operator 
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+ 5 t , t ,E%]XZ(t') + ^ dr ^(r)X£(r) j , (2) 

where H c0 — J2 a ec S a ^aX" a is the local cluster Hamil- 
tonian, T represents the imaginary time ordering opera- 
tor, and we used the notation (a/3) = /i. The hybridiza- 
tion and the effective magnetic field h% are the Weiss 
fields describing the effects of the rest of the system on 
the cluster. 

A cluster DMFT approach to a lattice problem consists 
of two elements: 1) a recipe for expressing the Weiss fields 
in terms of cluster quantities , i.e. a self-consistency con- 
dition, and 2) a recipe for determining lattice quantities 
from the relevant cluster counterparts, i.e. a periodiza- 
tion procedure. To carry out the first step we follow 
the CDMFT approach and construct a super-lattice by 
translating a cluster C so as to cover the original lattice 
and treat the cells as "single" sites with internal degrees 
of freedom. The self-consistency equation for A^ is de- 
termined by the condition that the Green's function for 
the Hubbard operators be the same for a cell and for the 
impurity cluster: 



J2 [m; 1 --e k ]- 1 = [m- i -a--e}-\ 
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where we used the tensor notation A = with a and 
b labeling the sites inside the cluster of size N c , and 
/i = (a,/3), v = (a',{3'). In Eq.J3J) M c is the irreducible 
two-point cumulant of the cluster defined as the sum of 
all two point diagrams generated by the strong coupling 
expansion of Eq. (J2J that are irreducible with respect to 
E and A, E represents the coupling constant matrix and 
the k summation is performed over the reduced Brillouin 
zone associated with a super-lattice with cells of size N c . 
The Weiss field h% is determined by 

E E cw) = k + EE A "ao) m, w 

i£C v beC v 

where the hybridization is evaluated at zero frequency. 
Within our approach, the cluster problem defined by 
Eq. J5J has to be solved self-consistently together with 
Eqs. 10) and (@}. Notice that only cluster quantities, 
in particular the irreducible cumulant M c , enter the self- 
consistency loop. 

The impurity model delivers cluster quantities and, to 
make connection with the original lattice problem, we 
need to infer from them estimates for the lattice Green's 
function. A natural way to produce these estimates is 
by considering the super-lattice construction described 
before and averaging the relevant quantities, which we 
denote below by W, to restore periodicity, namely 

W(i-j)*j^Yl W k,k+i-j> ( 5 ) 
s k 

where W and W SL are the lattice and super-lattice quan- 
tities, respectively, and N s represents the total number 
of sites. We stress that Eq. (|5j) represents a super-lattice 
average, not a cluster average. In particular, if W is the 
irreducible cumulant, all the contributions with k and 
k + i — j belonging to different cells are zero by con- 
struction. One possibility Q, is to periodize the Green's 
function 

G(k,c) = ± £ [M- 1 - EkU 1 e*(*.-*0 ( 6 ) 

C a,beC 

Ek being the Fourier transform of the "hopping" on the 
super-lattice, and N c the number of sites in a cell. A sec- 
ond possibility, suggested by the the strong coupling ap- 
proach investigated in this letter, is to first periodize the 
irreducible cumulant and then use it to reconstruct the 
lattice Greens function G(k,w) = [M _1 (k,w) - Ek] -1 . 
For example, within a four-site approximation (plaque- 
tte) we obtain after performing the average (JSJ and then 
taking the Fourier transform, 

M(k,w) = Mo(w)+Mi(a;) a(k) +M 2 H /3(k), (7) 

where a(k) = cos(k x ) +cos(fc y ), /3(k) = cos(k x ) cos(k y ) 
and Af p= {o,i,2} represents the on-site, nearest neighbor, 
and next nearest neighbor cluster cumulant, respectively. 



To test the dependence of the approach on the super- 
lattice construction, we also introduce an alternative self- 
consistency condition that involves the periodized lattice 
quantities, instead of the cluster quantities that appear 
in Eq. ©, in the spirit of PCDMFT0 but satisfying an 
explicit cavity construction. We define the hybridization 
function A as the sum of all the contribution to the clus- 
ter irreducible cumulant coming from outside the cluster 
and being connected to bare cumulants inside the cluster 
by two "hopping" lines, namely: 

(lUJ n ) = E aA K AB ( lWll )E S6 , (8) 

A,B 

where we used the matrix notation W a fc = and the 
matrix multiplication over \i and v is implied. In Eq. 
(JSJ, K^s(«w„) represents the cavity propagator, i.e. a 
Green's function which does not contain contributions 
arising from irreducible cumulants having at least one 
site index inside a certain cluster C of size N c . We as- 
sume that My has a finite range — rj\ < R, so that 
the terms that we subtract from the lattice cumulants to 
construct the cavity form a matrix M* which is nonzero 
only inside an extended cluster C ex t containing sites that 
can be coupled with the original cluster by a non-zero 
cumulant. Explicitly, M* 3 = My if at least one of the in- 
dexes belongs to the cluster C and zero otherwise, which 
insures that M*- is contained in C ex t- The propagator 
K.AB(iu n ) is a generalization of the cavity functional 
and we are interested to express it in terms of the lattice 
Green's function for sites (A, B) that can be connected 
with the cluster C via a hopping line. We assume that 
(A, B) E C ex t, i-e. the hopping has the same range as M y 
or smaller. For the extended cluster the cavity propaga- 
tor can be written as 

K = G - HM*[i + f M*] _1 H, (9) 

where all the matrices ~Wab are defined on the extended 
cluster, A,B <E C ex t, G is the lattice Green's function 
corresponding to the irreducible lattice cumulant M, 
Uab = {GM.~ 1 )ab, and £ AB = (EH -1 )ab. The lat- 
tice Green's function can be expressed directly in terms 
of cluster cumulants using Eq. (7J as 

G- X (k,w) = [Mo(w)+Mi(w) a(k)+M 2 H ^(k)]- 1 -^. 

(10) 

We note that equation JHJ together with Eq. © and 
(jlOJI can be viewed as an independent cluster scheme 
with variations that can be generated using different 
choices for the cavity matrix M*. This alternative clus- 
ter method, allows us to check that our results are self 
consistent by reintroducing the lattice cumulant into 
the DMFT equations. This is important since previ- 
ous results of a straightforward strong coupling expansion 
were shown to disappear in a more sophisticated DMFT 
treatment 0. 
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FIG. 1: Kinetic energy of the half-filled one-dimensional Hub- 
bard model as a function of the on-site interaction U at zero 
temperature using: the Bethe ansatz (red line), the cluster 
Green's function (black circles) and the lattice Green's func- 
tion obtained by periodizing G (green squares) , the self-energy 
(blue triangles) and the irreducible cumulant (magenta line). 
The results are obtained using two site CDMFT and an exact 
diagonalization (ED) impurity solver. 



We benchmark our approach, as in ref. 10], by com- 
puting the kinetic energy of the half filled one dimen- 
sional Hubbard model which is given known exactly from 
the Bethe ansatz. We also make a comparison with 
the alternative periodization procedures involving the 
Green's function and the self-energy. Shown in Fig. 
is the kinetic energy of the half-filled one-dimensional 
Hubbard model. The exact result from the Bethe ansatz 
(red) is used as a benchmark. We notice that the values 
of the kinetic energy given by the cluster Green's function 
(black) are significantly different from the exact result, 
while the curves obtained using the lattice Green's func- 
tion, extracted using various procedures, cluster around 
the Bethe ansatz line. We notice that the results ob- 
tained by periodizing the Green's function (green) and 
those obtained by periodizing the cumulant (magenta) 
are remarkably similar, especially in the strong coupling 
regime. We observed a very similar behavior in the two 
dimensional case. We conclude that within CDMFT and 
other related cluster schemes, when applied to small clus- 
ters, observables should be always extracted from the 
physical lattice quantities and not from their cluster coun- 
terparts. Because in our generalized strong-coupling con- 
struction of the cluster approximations hopping is treated 
on equal footing with other non-local contributions to the 
Hamiltonian, such as spin-spin interaction, the conclu- 
sions derived from the calculation of the kinetic energy 
extend to all non-local physical quantities, for example 
to the spin-spin correlation function. In contrast with 
non-local quantities, the local physical quantities are well 
approximated by their cluster values which have to be 
preserved by the reconstruction schemes. For a homo- 
geneous cluster (as, for example, the link or the plaque- 



tte), periodizing the Green's function automatically sat- 
isfies this condition for all one-particle quantities as, by 
construction, G\f tt = G c aa . The cumulant periodization 
scheme also generates a local Green's function in good 
agreement with G c . However, the self-energy scheme fails 
at half filling and for small doping values as it generates 
spurious states in the gap. 

As a first application of our method to a strongly cor- 
related metal, we study the two-dimensional Hubbard 
model using a four-site cluster approximation. In gen- 
eral, the lattice Green's function can be written as 



G(k,u/) = 



1 



lu — r(k, lu) — ir)(k, to) ' 



(11) 



where lu) represents the imaginary part of the self- 
energy and r(k, lu) = e(k) — \i + i?eE(k, lu) is the energy. 
In the self-energy periodization scheme doping values, 
E(k, lu) is a linear combination of the lattice self-energies 
given by 

S(k, lu) = E (w) + Ei(w) a(k) + E a (w) /3(k). (12) 

In the cumulant re-construction scheme, which describes 
better the system near the Mott transition, the lattice 
self-energy is given by a highly non-linear relation 



E(k,w) = 



LU — /! 

LU + fl - Y.A 



(13) 



t (l-a + P) \{l + a + (3) ~ 
x) + ^-E_b w + p - Sc . 



where a(k) and (3(k) were defined above, and the diago- 
nal cluster self-energies are = E — E 2 and Tibic) — 
E T 2Ei + E 2 . Using an exact diagonalization as a 
CDMFT impurity solver^l] one finds that at zero tem- 
perature the imaginary parts of the cluster self-energies 
go to zero at zero frequency. For the real parts, on the 
other hand, we distinguish two regimes. At large dop- 
ings the diagonal cluster self-energies are dominated by 
the local component Eo and Eq. I|14|l reduces in the first 
approximation to Eq. (|12[) . In this regime the physics 
is almost local with small corrections due to short-range 
correlations. All the periodization schemes converge and 
the single-site DMFT represents a good first order ap- 
proximation. In contrast, close to the Mott transition 
the short-range correlations become important and the 
off-diagonal components of the cluster self-energy become 
comparable with Eo. As a consequence, at zero frequency 
the denominators in Eq. (|14|l may acquire opposite signs 
generating a divergence in the lattice self-energy. This 
pole of E(k, lu = 0), or equivalently of r(k), gives rise to 
a zero of the lattice Green's function. We show in Fig. [21 
the renormalized energy, r(k), and the spectral function 
A(k, u> = 0) = — l/7r7roG(k, 0) for a two-dimensional 
Hubbard model with U = 8t at zero temperature for two 
values of doping. For n — 0.78 (left panels) we have a 
large electron- type Fermi surface (blue line in the r(k) 
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FIG. 2: Renormalized energy, r(k), (upper panels) and spec- 
tral function, A(k), (lower panels) for the 2D Hubbard model 
with U = 8t and T = 0. The color code for the upper panels 
is: green (r < 0), blue (r = 0), yellow (r > 0), red (r — * oo). 



panel) separating the occupied region of the Brillouin 
zone (green), defined by r(k) < from the unoccupied 
region (yellow) defined by r(k) > 0. The Fermi surface 
can be also traced in the A(k) panel as the maximum of 
the spectral function. On the other hand, for n = 0.92 
a qualitatively different picture emerges. The Fermi sur- 
face (blue line) is now represented by a hole pocket and, 
in addition, we have a line of zeros of the Green's function 
(red line) close to the (it, ir) region of the Brillouin zone. 
Furthermore, there is no one-to-one correspondence be- 
tween the Fermi surface and the maximum of the spectral 
function. This behavior has two origins: 1) the proximity 
of a zero line suppresses the weight of the quasiparticle 
on the far side of the pocket, and 2) for k-points corre- 
sponding to r(k) ^ the quasiparticles are pushed away 
from to = and a pseudogap opens at the Fermi level. 
We show this explicitly in Fig. 01 by comparing the low 
frequency dependence of the spectral function in three 
different points of the Brillouin zone, marked by A, B and 
C in Fig. [21 Notice the suppression of the zero frequency 
peak at point B and the frequency shift S = — 0.05t of 
the peak at point C. The cumulant approach provides 
a simple interpretation of this effect, observed in photo 
emission experiments m terms of the emergence of 
infinite sclf-cncrgy lines or cquivalcntly Luttingcr lines ( 
lines of zeros of the Greens function). 

In conclusion, our strong coupling CDMFT study of 
the Hubbard model shows that the lightly doped system 
is characterized by a small, closed Fermi line that appears 
in the zero frequency spectral function as an arc due to 
the presence of a line of zeros of the Green's function near 
the "dark side" of the Fermi surface. These lines appear 
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FIG. 3: Frequency dependence of the spectral function for 
three points in the Brillouin zone marked by A, B, and C in 
Fig. 



near the Mott insulator in order to satisfy the generalized 
Luttinger theorem[l3|. 

N f d 2 k 

— = 2 - — (14) 

V JReG(k )W =o)>o ( 27r ) 

While this theorem is not exactly satisfied in our finite 
cluster calculation due to the small cluster approxima- 
tion, even this small cluster size has all the qualitative el- 
ements needed to interpret the physics of the underdoped 
regime which is controlled by the position of the lines of 
zeros and infinities. Notice that the standard Luttingcr 
theorem will be always violated in a system character- 
ized by a pseudogap, and it is strongly violated in our so- 
lution as the area contained in the pockets shrinks. The 
tendency to increase the area contained by the Luttinger 
surface as the pockets shrink is correctly captured by our 
approximation. Furthermore, we identify the pseudogap 
which is seen in leading edge study of photo emission 
experiments as the small negative shift of the spectral 
weight in points of the Brillouin zone that are not on the 
Fermi line (for example point C in Fig. |5J), which is dis- 
tinct from the larger gap between the peaks above and 
below the Fermi level (see Fig. EJ. Remarkably, the lines 
of poles of the self energy appear first far from the Fermi 
surface. This is a strong coupling instability which has 
no weak coupling precursors on the Fermi surface. Our 
results raise an interesting question. If the evolution in 
Figure [3 from large to small doping is continuous, it has 
to go through a critical point where the topology of the 
Fermi surface (and perhaps that of the lines where the 
self energy is infinite) changes. This topological change 
and its possible connection to an underlying critical point 
at finite doping in the cuprate phase diagram deserves 
further investigation. 
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